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STATEMENT  OF  THE  PROBLEM. 


Epoxy  resins  containing  the  bisphenoI-A  unit  are  among  the  most 
commonly  employed  crosslinked  or  network  polymer  systems.  These 
materials  display  significant  relaxatons  below  the  glass  transition.  These 
sub-glass  transition  relaxations  of  bisphenoI-A  based  epoxide  resins  have 
been  widely  studied  by  dynamic  mechanical  responsei-i6,  dielectric 
responsei7-i8  and  by  solid  state  NMR19*22.  in  both  the  mechanical  and 
dielectric  studies,  the  principal  low  temperature  loss  peak  (  the  p  peak)  Is 
quite  broad  and  lies  between  -80  and  -40  ©C  at  a  frequency  of  1  Hz.  The 
effects  of  crosslink  density  and  changes  In  chemical  structure  on  the  p 
peak  have  been  studied. 

The  interpretation  of  the  main  loss  peak  in  terms  of  a  molecular 
level  motion  remains  unclear.  One  group  of  investigators6*7,23  has 
attributed  this  loss  to  motion  of  the  diphenylpropane  unit,  specifically 
phenylene  group  rotationis.  Others  attribute  the  relaxation  to  motion  of 
the  hydroxyether  groups3-4, 8-12.22.24  while  some  of  the  more  recent 
mechanical  studies  feel  the  p  peak  is  a  composite  of  motions  of  both  of 
these  unitss.  13*17.  According  to  this  latter  interpretation,  the  low 
temperature  side  (-110  to  -8O0C)  of  the  p  peak  arises  from  motion  of  the 
phenylene  groups  in  the  diphenyl  propane  and  the  high  temperature  side  of 
the  p  peak  arises  from  the  motion  of  the  hydroxyether  groups  (-  70  to  -55 
oC).  To  the  current  authors,  the  literature  is  not  definitive  on  the 
molecular  sources  of  the  sub-glass  transition  motion.  It  is  difficult  to 
establish  the  molecular  origin  of  a  mechanical  loss  by  observing  changes 
In  mechanical  loss  produced  by  structural  variations.  It  is  the  purpose  of 
this  contribution  to  use  solid  state  deuterium  echo  line  shapes  to  more 
clearly  identify  the  structural  source,  geometry,  time  scale  and 
temperature  dependence  of  the  molecular  motions  present  and  to 
associate  them  with  the  mechanical  response  of  the  bisphenol  epoxide 
resin. 

Deuterium25-27  and  carbon-13  28  |jne  shape  studies  have  been 
established  as  a  definitive  technique  for  the  characterization  of  motion  in 
solid  polymers  and  for  establishing  the  relationship  between  molecular 
motions  and  the  bulk  mechanical  response  of  the  solid  material.  To  apply 
these  techniques  most  directly,  isotopically  labeled  samples  must  be 
prepared.  For  the  current  study,  the  obvious  sites  to  label  with  deuterium 
in  separate  samples  are  the  phenylene  rings  of  the  diphenylpropane  portion 
and  the  methylene  groups  of  the  hydroxyether  portion.  This  will  allow  for 


detailed,  direct  monitoring  of  motion  In  each  of  the  two  key  regions  of  the 
material  separately.  To  test  the  effects  of  crosslinking  or  curing  the 
dynamics  of  a  cured,  crosslinked  material  will  be  compared  with  a  linear 
material. 


SUMMARY  OF  RESULTS. 


Experimental 

Four  deuterium  labeled  materials  were  prepared  for  this  study.  First 
linear  poIy(diglycidyl  ether  of  bisphenol  A)  (poly(DGEBA))  was  prepared 
with  deuterated  methylene  units  in  the  hydroxyether  portion.  A  second 
sample  of  the  same  linear  polymer  was  prepared  with  deuterated 
phenylene  units  in  the  diphenylpropane  portion.  The  third  labeled  sample 
is  made  by  curing  the  methylene  labeled  DGEBA  with  4,4’  methylene 
dianlline  (MDA)  to  produce  a  model  one  to  one  network.  A  similar  network 
sample  is  prepared  by  curing  the  phenylene  labeled  DGEBA  with  MDA.  The 
crosslinked  resins  were  synthesized  by  the  addition  of  1  equivalent  of 
MDA  to  4  equivalents  of  DGEBA.  They  were  then  cured  by  heating  for  16 
hours  at  55oC  followed  by  2  hours  at  125oC  and  2  hours  at175oC.  The  linear 
polymers  were  made  by  addition  of  BIsphenol-A  and  DGEBA  (1:1)  followed 
by  the  same  cure  conditions.  The  structures  of  the  four  labeled  systems 
are  shown  in  Figure  1. 

The  deuterium  NMR  experiments  were  carried  out  with  a  Bruker  MSL- 
300  spectrometer  operating  at  a  deuterium  resonance  frequency  of  46.07 
MHz.  Temperature  was  regulated  with  a  Bruker  VT-1000  control  unit.  The 
deuterium  spectra  were  obtained  using  a  quadrupole  echo  pulse  sequence 
with  echo  delay  times  of  15,  30,  50,  80,  and  120  \is.  The  90°  pulse  width 
was  2.5p.s  and  the  recycle  time  of  the  experiment  was  5  to  20  s  depending 
on  the  sample  and  the  temperature. 

Results 

Solid  echo  spectra  on  the  four  samples  as  a  function  of  echo  time 
and  temperature  over  the  range  of  about  -100  to  lOOoC  are  shown  In 
Figures  2,  3,  4,  and  5.  At  the  lowest  temperatures  (-80oC),  Pake  patterns 
are  observed  corresponding  to  nearly  rigid  specimens.  The  quadrupole 
coupling  constant  for  the  methylene  labeled  systems  is  170  kHz;  and  for 
the  phenylene  labeled  systems,  180  kHz.  As  temperature  is  raised,  the 
methylene  labeled  systems  partially  collapse  to  a  line  shape 
corresponding  to  an  asymmetric  tensor  which  has  generally  been 
associated  with  the  onset  of  trans-gauche  Isomerization  plus  rapid 
libratlon26-27.  The  phenylene  labeled  systems  partially  collapse  as 
temperature  Is  raised  and  show  a  splitting  at  high  temperatures  which  is 
indicative  of  n  flips  plus  rapid  librationss.  Both  methylene  and  phenylene 
labeled  systems  show  lineshape  changes  as  a  function  of  echo  delay  time 
at  certain  temperatures  which  verifies  the  presence  of  motions  on  time 
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scales  comparable  to  the  selection  of  echo-  delay  times.  Also  the  spectra 
show  significant  intensity  reductions  with  increasing  echo  delay  time  at 
these  same  temperatures  further  identifying  the  time  scale  of  motion  as 
comparable  to  the  echo  delay  timesss. 

Interpretation 

To  quantitatively  determine  the  rates  and  amplitudes  of  motion,  the 
deuterium  line  shapes  are  simulated  according  to  the  procedures  described 
by  Wittebort29  and  Mehringso.  Simulation  calculations  were  performed  on  a 
DEC  UNIX  workstation  and  the  simulated  spectra  Include  the  effects  of 
finite  pulse  width.  Because  of  the  complexities  of  molecular  dynamics  in 
polymeric  glasses  several  additions  to  these  procedures  must  be  made. 

The  first  Is  that  motion  In  a  glass  cannot  be  characterized  by  a  single  rate 
but  rather  the  motion  Is  heterogeneous  and  must  be  described  by  a 
distribution  of  ratesi9i3i.  in  this  regard  one  of  the  more  successful 
functions  for  describing  relaxation  In  polymers  is  the  stretched 
exponential  or  Kohlrausch-  William-  Watts  functloni 9.32-32 


(})(t)  =  exp  -(t/Tp)“  (1) 

which  can  be  written  as  sum  of  exponentlals33.  The  fractional  exponent,  a, 
controls  the  breadth  of  the  distribution  and  the  characteristic  time,  tp, 
controls  the  time  scale.  Accordingly  these  two  parameters  will  be  quoted 
when  summarizing  the  rates  or  time  scales  of  motion. 

Two  jump  motions  are  considered  In  the  interpretation:  trans-gauche 
isomerization  for  the  methylene  units  and  k  flips  for  the  phenylene  units. 
However  In  disordered  glasses  a  distribution  jump  angles  is  required  for 
both  of  these  motions  and  following  EngIlsh26-27  ©t  ©i  ©  Gaussian 
distribution  of  jump  angle  centered  around  the  conformational  minimum 
will  be  used.  The  usual  Gaussian  breadth  parameter  a  is  reported  to 
characterize  this  aspect  of  the  dynamics. 

The  last  dynamic  influence  on  the  line  shape  spectra  comes  from  the 
presence  of  llbratlonal  motion  which  Is  relatively  rapid  on  the  line  shape 
collapse  time  scale.  To  include  the  effect  of  this  motion,  the  quadrupole 
tensor  is  pre-averaged  before  the  line  shape  collapse  calculation  is 
mad 625-28.  Again  a  distribution  of  llbratlonal  amplitudes  Is  often  required 
in  gIasses26-27  but  to  simplify  the  calculations  and  reduce  the  time 
required  on  the  computer  a  single  librational  amplitude  was  employed. 
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Improved  fits  could  be  achieved  if  a  distribution  of  librational  amplitudes 
was  employed  but  the  essential  aspects  of  the  dynamics  are  captured 
without  a  distribution  for  this  influence  on  the  line  shape.  The  amplitude 
of  libration  is  temperature  dependent  in  the  modeling  presented  here  and 
is  allowed  to  increase  with  temperature.  Improved  simulations  could  also 
be  produced  if  the  librational  motion  was  given  an  actual  rate  and  not 
simply  assumed  to  be  fast  on  the  NMR  time  scale.  Again  this 
simplification  is  made  to  limit  the  simulation  to  reasonable  times  on  a 
work  station. 

To  properly  calculate  the  changes  in  amplitude  with  increasing  echo 
delay  times,  the  spin-spin  relaxation  time,  T2,  must  be  included  and  It  can 
be  estimated  from  the  low  temperature  spectra  or  the  high  temperature 
spectra  where  the  reduction  of  intensity  associated  with  dynamics  are 
minimal. 

The  calculated  line  shapes  and  reductions  factors  are  compared  with 
the  experimental  results  in  Figure  2  for  the  linear  poly(DGEBA)  with 
methylene  deuterated  hydroxyether  units.  The  simulation  parameters  are 
summarized  in  Table  1.  The  analogous  comparisons  of  experiment  and 
simulation  for  the  crosslinked  methylene  deuterated  hydroxyether  system 
are  show  in  Figure  3  and  the  associated  simulation  patterns  are  given  In 
Table  2.  Figure  4  contains  the  comparison  of  experiment  and  simulation 
for  the  phenylene  labeled  linear  poly(DGEBA)  and  the  simulation 
parameters  are  contained  in  Table  3.  The  last  set  of  simulations  for  the 
phenylene  labeled  cross  linked  system  are  compared  with  experiment  in 
Figure  5  and  the  parameters  are  listed  in  Table  4. 

Apparent  activation  energies  are  calculated  from  a  fit  of  the 
simulation  values  of  tp  according  to  the  Arrhenius  equation 

Xp  =  x«  exp(Ea/RT)  (2) 

which  are  shown  In  Figure  6.  Since  the  time  scales  of  motion  are  all 
rather  similar,  a  fixed  value  of  the  prefactor  of  10-i4  is  used  in  all  fits. 
Table  5  gives  a  summary  of  the  apparent  activation  energies  obtained  for 
the  motion  in  each  system. 

The  amplitude  of  libration  as  a  function  of  temperature  is  shown  in 
Figure  7.  Actually  the  libration  amplitude  has  a  nearly  linear  dependence 
on  the  square  root  of  temperature  as  has  been  observed  In  other  systems 
so  Figure  7  is  plotted  accordingly. 


Discussion 


The  most  rapid  motion  observed  in  the  four  systems  studied  is  for  tc 
flips  of  the  phenylene  group  in  the  linear  poly(DGEBA)  and  the  second  most 
rapid  motion  is  for  tc  flips  of  the  phenylene  groups  In  the  crosslinked 
network.  The  effect  of  crosslinking  on  phenylene  group  motion  is  about 
equivalent  to  a  shift  of  time  scale  corresponding  to  a  temperature  change 
of  5  or  10  degrees.  The  breadth  parameter  determining  the  distribution  of 
correlation  times,  a,  Is  0.6  for  the  linear  system  and  0.7  for  the 
crosslinked  network.  This  corresponds  to  a  fairly  narrow  distribution  of 
correlation  times  and  a  prediction  of  the  temperature  of  the  maximum  of 
the  mechanical  loss  peak  can  be  obtained  by  extrapolating  the  fit  of  the 
NMR  data  to  a  time  scale  of  1  Hz.  The  1Hz  line  is  indicated  In  Figure  6  and 
the  corresponding  temperature  Is  -95oC  for  the  linear  resin  as  listed  in 
Table  5.  The  position  of  the  mechanical  loss  peak  predicted  from  the  NMR 
data  for  the  crosslinked  network  is  11  degrees  higher.  The  positions  of 
these  predicted  mechanical  loss  peaks  place  them  on  the  low  temperature 
side  of  the  p  peak  but  part  of  the  overall  envelope  of  loss  associated  with 
P  peak.  This  is  consistent  with  the  views  of  some  of  the  mechanical 
investlgators8.i3*i7  though  mechanical  data  Is  so  featureless  that  changes 
produced  by  structure  variation  leave  considerable  uncertainty  as  to  the 
position  of  the  peak  to  be  associated  with  phenylene  group  motion.  The 
absolute  uncertainty  of  the  NMR  estimate  of  the  temperature  of  the  loss 
peak  is  between  5  and  10  degrees  though  the  error  In  the  relative 
difference  of  1 1  degrees  between  linear  and  crosslinked  Is  closer  to  2 
degrees. 

The  activation  energy  of  48-51  kJ/mole  for  the  %  flip  motion  is 
typical  of  those  found  for  phenylene  groups  in  the  backbones  of  polymers. 
For  instance,  the  apparent  activation  energy  for  the  same  motion  in  linear 
polycarbonate  is  40  to  50  kJ/mole28,3i.  in  Garroway’s  NMR  studyis,  a  value 
of  60kJ/mole  is  obtained  for  phenylene  group  motion  in  a  crosslinked 
epoxy. 

The  trans-gauche  Isomerization  of  the  methylene  groups  in  the 
hydroxyether  portion  is  a  slower  motion  than  tc  flips  in  either  the 
crosslinked  or  linear  systems  though  the  two  geometries  of  motion  are 
only  separated  by  about  a  decade  in  time  or  frequency.  The  effect  of 
crosslinking  is  to  slow  the  isomerization  motion  slightly  relative  to  the 
linear  system.  The  slowing  corresponds  to  a  temperature  shift  of  about  5 
degrees.  The  apparent  activation  energy  obtained  for  this  motion  is  55-58 
kJ/mole  which  leads  to  a  slightly  larger  temperature  separation  at  the 
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time  scale  of  1  Hz  corresponding  to  mechanical  measurements.  Since  the 
distributions  of  correlation  times  are  again  relatively  narrow  the 
relaxation  map  of  Figure  6  leads  to  a  predicted  temperature  for 
mechanical  loss  of  -70  oc  for  the  linear  system  and  -57oC  for  the 
crosslinked  system.  This  difference  between  linear  systems  and 
crosslinked  systems  seems  reasonable  relative  to  mechanical  data  on 
similar  systemsie.  These  temperatures  fall  on  the  high  side  ofjhe 
envelope  of  mechanical  loss  associated  with  the  p  peak.  Thus  in  both  the 
linear  and  crosslinked  systems  the  two  motions,  %  flips  and  trans-gauche 
isomerization,  would  likely  overlap  In  a  mechanical  experiment  and 
contribute  to  the  overall  mechanical  peak  described  as  the  p  peak.  This 
coincidence  between  the  mechanical  experiment  and  the  prediction  based 
on  the  interpretation  of  the  NMR  lineshape  data  can  be  further  visualized 
when  a  calculatedsi  loss  peak  is  plotted  and  compared  with  the 
experimental  loss  data22  (Figure  8).  The  equationsi  used  to  calculate  the 
shear  loss  from  the  rate  and  temperature  dependence  established  from  the 
interpretation  of  the  NMR  data  is 


G(<or  =  isin(o>t)[-f(t)]dt 

0 


(3) 


The  calculation  provides  information  on  the  breadth  of  the  loss  peak  and 
the  position  in  frequency  and  temperature  but  not  on  the  amplitude  of  the 
loss  peak.  An  amplitude  Is  assumed  which  matches  experimental 
observation.  The  calculated  loss  peak  Includes  the  low  temperature  (fast) 
ring  motion  and  the  slower,  higher  temperature  chain  motion  and  these  can 
be  discerned  in  the  calculated  line.  It  is  clear  from  the  comparison  that 
the  p  relaxation  contains  contributions  from  further  slower  motions  not 
detected  in  our  NMR  study  which  contribute  to  loss  at  the  high 
temperature  side  of  the  mechanical  relaxation.  There  are  two  potential 
sources  for  this  slower  motion  in  the  crosslinked  system.  One  is  motion 
associated  with  the  crosslinking  agent  Itself,  MDA.  Since  a  sample  with 
deuterium  labeled  MDA  is  not  available  we  would  not  observe  any  local 
motion  of  the  MDA  moiety  in  the  resin.  The  crosslinked  epoxy  might  also 
exhibit  an  additional  slow  motion  component  associated  with  the  presence 
of  the  crosslinker  at  one  side  of  the  methylene  units  relative  to  the  more 
symmetric  linear  system.  This  could  explain  the  slightly  lower  quality  of 
the  simulations  for  the  crosslinked  system  (figures  3  and  5).  The  fit  could 


obviously  be  improved  by  the  inclusion  of  more  parameters  to  account  for 
this  but  this  does  not  seem  justified  in  view  of  the  inherent  information 
content  of  the  iineshapes  and  smce  the  general  interpretation  is  not 
enhanced. 

The  apparent  activation  energy  of  ca.  57  kj/mole  is  high  relative  to 
other  observed  trans-gauche  isomerization  energies  for  methylene  units. 
For  instance  the  value  obtained  by  English  et  al  on  the  methylene  groups  in 
nylon  Is  33kJ/mole.  The  higher  value  may  reflect  the  presence  of  the 
hydroxyl  groups  on  this  portion  of  the  network  which  can  participate  in 
hydrogen  bonding. 

Conclusion 

Deuterium  line  shape  data  and  the  subsequent  analysis  identifies  the 
role  of  K  flips  of  the  bisphenol  A  unit  and  trans-gauche  isomerization  of 
the  mehtylene  units  in  the  mechanical  response  of  epoxy  networks.  Both 
motions  contribute  to  the  main  sub-glass  transition  relaxation  process, 
the  p  process,  with  the  7t  flips  occuring  at  lower  temperatures  and  the 
trans-gauche  isomerization  at  higher  temperatures  but  with  overlap  of 
the  two  relaxations.  The  activation  energy  of  the  trans-gauche 
isomerization  Is  about  7kJ/mol  higher  than  the  n  flip  process.  The  effect 
of  cross-linking  is  a  shift  of  either  relaxation  process  about  ten  degrees 
at  the  typical  frequency  of  a  mechanical  experiment  of  1  Hz.  Cross-linking 
may  produce  some  subtle  changes  in  the  distribution  of  relaxation  times 
but  It  was  not  necessary  to  include  these  to  produce  reasonable 
simulations. 
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TABLE  I 

Simulation  parameters  for  the  chain  (methylene) 
motion  for  the  linear  epoxy. 


T(oC) 

xp(s) 

a 

Libratlon(o) 

0(0) 

T2(ms) 

70 

6,6x10-6 

0.6 

85 

18 

1.5 

48 

1.6x10-5 

0.6 

80 

18 

,1.5 

35 

2.5x10-5 

0.6 

76 

18 

1.5 

24 

4.5x10-5 

0.6 

71 

18 

1.5 

0 

1.5x10-4 

0.6 

61 

18 

1.5 

-20 

5.0x10-4 

0.6 

52 

18 

1.5 

TABLE  II 

Simulation  parameters  for  the  chain  (methylene)  motion 
for  the  crossllnked  epoxy. 


T(oC) 

Tp(S) 

a 

Libratlon(o) 

o(o) 

T2(ms) 

90 

1.0x10-5 

0.7 

75 

18 

0.8 

70 

1.5x10-5 

0.7 

72 

18 

0.8 

48 

3.5x10-5 

0.7 

65 

18 

0.8 

35 

6.0x10-5 

0.7 

60 

18 

0.8 

24 

1.4x10-4 

0.7 

56 

18 

0.8 

0 

6.0x10-4 

0.7 

47 

18 

0.8 

-20 

2.8x10-3 

0.7 

39 

18 

0.8 
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TABLE  III 

Simulation  parameters  for  the  ring  motion  for  the  linear  epoxy. 


T(oC) 

Tp(s) 

a 

Libration(o) 

a{o) 

T2(ms) 

70 

2.4x10-7 

0.6 

50 

25 

0.28 

48 

6.5x10-7 

0.6 

47 

25 

0.28 

35 

1.4x10-6 

0.6 

44 

25 

0.28 

24 

2.8x10-6 

0.6 

41 

25 

0.28 

0 

1.5x10-5 

0.6 

35 

25 

0.28 

-20 

4.0x10-5 

0.6 

30 

25 

0.28 

Simulation 

parameters 

TABLE  IV 

!  for  the  ring  motion  for  the 

crossiinked 

epoxy. 

T(oC) 

xp(s) 

a 

Libration(o) 

a(o) 

T2(ms) 

90 

2.9x10-7 

0.7 

50 

25 

0.20 

70 

6.7x10-7 

0.7 

47 

25 

0.20 

48 

1.8x10-6 

0.7 

43 

25 

0.20 

35 

4.1x10-6 

0.7 

40 

25 

0.20 

24 

7.8x10-6 

0.7 

37 

25 

0.20 

0 

3.7x10-5 

0.7 

31 

25 

0.20 

-20 

2.8x10-4 

0.7 

26 

25 

0.20 

Activation  parameters 

TABLE  V 
for  linear  and 

crossiinked  epoxies 

Ea(kJ/moi) 

Temp  at  IHz(oC) 

Linear,  chain  motion 

54.6 

-70 

Crossiinked,  chain  motion 

57.9 

-57 

Linear,  ring  motion 

47.7 

-95 

Crossiinked,  ring  motion 

50.7 

-84 
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Figure  captions. 


Figure  1: 
Figure  2: 

Figure  3: 

Figure  4: 

Figure  5: 

Figure  6: 

Figure  7: 
Figure  8: 


Chemical  structures  of  the  four  deuterium  labeled  epoxy 
systems. 

2H  experimental  (left)  and  simulated  (right)  lineshapes  for  the 
chain  labeled  linear  epoxy  as  a  function  of  echo  delay  times: 

15,  30,  50,  80,  120ps  from  bottom  to  top;  reduction  factors 
are  listed  at  the  right. 

2H  experimental  (left)  and  simulated  (right)  lineshapes  for  the 
chain  labeled  crosslinked  epoxy  as  a  function  of  echo  delay 
times:  15,  30,  50,  80,  120ps  from  bottom  to  top;  reduction 
factors  are  listed  at  the  right. 

2H  experimental  (left)  and  simulated  (right)  lineshapes  for  the 
ring  labeled  linear  epoxy  as  a  function  of  echo  delay  times: 

15,  30,  50,  80,  120ps  from  bottom  to  top;  reduction  factors 
are  listed  at  the  right. 

2H  experimental  (left)  and  simulated  (right)  lineshapes  for  the 
ring  labeled  crosslinked  epoxy  as  a  function  of  echo  delay 
times:  15,  30,  50,  80,  120)xs  from  bottom  to  top;  reduction 
factors  are  listed  at  the  right. 

Arrhenius  plot  for  both  the  chain  and  ring  motion  In  the  epoxy 
systems  studied.  The  rate  for  comparison  with  the  mechanical 
relaxation  (IHz)  is  shown  as  a  dotted  line. 

The  temperature  dependence  of  the  amplitude  of  the  libratlonal 
motion  of  the  chain  methylenes  and  the  phenylene  rings. 

A  comparison  of  the  experimental  mechanical  response22 
(circles)  with  the  prediction  (line)  based  on  the  ring  and  chain 
dynamics  obtained  from  the  NMR  analysis. 
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